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The Lieb-Mattis theorem about antiferromagnetic ordering of energy levels on bipartite lattices 
is generalized to finite-size two-leg spin-1/2 ladder model frustrated by diagonal interactions. For 
reflection-symmetric model with site-dependent interactions we prove exactly that the lowest en- 
ergies in sectors with fixed total spin and reflection quantum numbers are monotone increasing 
functions of total spin. The nondegeneracy of most levels is proved also. We also establish the 
uniqueness and obtain the spin value of the lowest-level multiplet in the whole sector formed by 
reflection-symmetric (antisymmetric) states. For a wide range of coupling constants, we prove that 
the ground state is a unique spin singlet. For other values of couplings, it may be also a unique 
^ I spin triplet or may consist of both multiplets. Similar results have been obtained for the ladder 

$_( , with arbitrary boundary impurity spin. Some partial results have also been obtained in the case of 

' periodical boundary conditions. 

^ I PACS numbers: 75.10.Jm, 75.40.Cx, 75.45. +j, 75.10.Pq, 75.50.Ee 
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^ , I. INTRODUCTION 

a , 

Nowadays, the frustrated spin systems are the subject of intensive study<ii2.i2.ii The interest on them is stimulated 
by recent progress in synthesizing corresponding magnetic materials*^ In these models, due to competing interactions, 
the classical ground state cannot be minimized locally and usually possesses a large degeneracy. The frustration 
can be caused by the geometry of the spin lattice or by the presence of both ferromagnetic and antiferromagnetic 
interactions. 

In contrast, the ground states of classical models on bipartite lattices are Neel ordered and are unique up to the 
global spin rotations. The bipartitness means that the lattice can be divided into two sublattices A and B, so that 
fS| ■ all interactions within the same sublattice are ferromagnetic while the interactions between different sublattices are 

■ antiferromagnetic. For the Neel state, all spins of A have same direction while all spins of B are aligned in opposite 
direction. 

The quantum fluctuations destroy Neel state and the quantum model, in general, has more complicated ground state. 
However, for bipartite spin systems, the quantum ground state inherits some properties of its classical counterpart. 
In particular, Lieb and Mattis proved that the quantum ground state of finite-size system is a unique multiplet with 
total spin Sgs = \Sa — Sb\, which coincides with the spin of the classical Neel state. Here, Sa and Sb are the highest 
• possible spins on corresponding sublattices.-'^i2^ This feature of bipartite systems looks natural, since in the limit of 
P5 I large spin values the quantum model approaches to the classical one. Moreover, a simple general rule that the energy 
O ■ increases with increasing spin, which is conventional for classical antiferromagnets, makes a sense in this case, too. 

More precisely, the lowest energy Es in the sector, where the total spin is equal to S, is a monotone increasing function 
^ , of the spin for any S > Sgs-^ This property is known as Lieb-Mattis theorem about antiferromagnetic ordering of 
r*"! ■ energy levels. Lieb-Mattis theorem is very important because it provides information about the ground state and 
rS ] spectrum of bipartite spin systems without exact solution or numerical simulation. Recently, it has been generalized 

■ to SU(n) symmetric quantum chain with defining representation and nearest-neighbor interactionsi^ A ferromagnetic 
■ ■ ■ ' ordering of energy levels has been formulated and proven also for XXX Heisenberg chains of any spin.-s^ 

More limited number of exact results is available for frustrated quantum spin systems. It is difficult to treat 
even their classical ground states as was mentioned above. The usual Lieb-Mattis theorem is no longer valid for 
frustrated systems. Recently, however, Lieb and Schupp proved rigorously that a reflection-symmetric spin system 
with antiferromagnetic crossing bonds possesses at least one spin-singlet ground state i"^"'^^ Moreover, under certain 
additional conditions, all ground states become singlets. A lot of frustrated spin systems satisfy the conditions of Lieb- 
Schupp approach. Two-dimensional pyrochlore antiferromagnet having checkerboard lattice structure is an example 
of such type of system. Another example is well-known two-dimensional antiferromagnetic Ji — J2 model, namely, 
the Heisenberg model on square lattice with diagonal interactions, which has been studied recently J^i^^ii^ii^ Note 
that this approach does not provide any information about the degeneracy of ground states. Thus, for frustrated spin 
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systems, this question still remains open. We mention also that the method used in Refs. [13 and [Tl| is restricted to 
the systems, which do not have a spin on their symmetry axis. 

At the same time, recently, certain signs in favor of Lieb-Mattis theorem for frustrated spin systems have appeared. 
For Ji — J2 model, the number of arguments based on exact diagonalization and spin-wave approximation had been 
presented in support of the fact that the antiferromagnetic ordering of energy levels is preserved under weak frustration 
caused by diagonal couplingsi^ Recentl y, f or the same model, the condition Sgs = \Sa — Sb\ for the ground state 
spin has been tested numerically in Ref. 'ITI. The authors came to the conclusion that it remains true provided that 
the frustration is sufficiently weak in order to destroy ground-state Neel order. The Lieb-Mattis property had been 
observed also at finite size spectrum of the Heisenberg model on the triangular and Kagome latticeSf^^ All these 
investigations suggest that in many cases, even a relatively strong frustration cannot destroy the antiferromagnetic 
ordering of energy levels. 

On the other hand, recently, the spin ladder systems have attracted much attention. The interest to them has 
been generated by significant progress made within the last years in fabrication of such type of compounds having a 
structure similar to the structure of two-dimensional high-temperature superconductors.— Ladder systems are simpler 
for study than their more complex two-dimensional counterparts. Moreover, the powerful theoretical and numerical 
methods, elaborated for one-dimensional models, can be applied to study them. The two-leg ladder frustrated by 
diagonal interactions is one of the simplest frustrated spin systems. It can be viewed as a quasi-one-dimensional 
analog of Ji — J2 model. Together with various generalizations, it has been investigated intensively during the 
last decadei ^°i^^'^^i^^'^"^i^^'^^'^^ The ground-state phase diagram consists of two topologically distinct Haldanc phases 
separated by the curve (or surface) of phase transition. 

Inspired by aforementioned activities, in this paper we generalize the Lieb-Mattis ordering theorem to two-leg spin- 
1/2 ladder model with diagonal interactions. We consider the system, which rests invariant under the reflection with 
respect to the longitudinal symmetry axis of the ladder. Up to our knowledge, this is the first exact result related to 
the energy level ordering of a frustrated spin system. 

In Sec. nil we formulate and prove an analog of Lieb-Mattis theorem for frustrated spin-1/2 ladder with free bound- 
aries. The reflection symmetry spits the total space of states into two invariant sectors, composed correspondingly 
from reflection-symmetric and reflection-antisymmetric states. We establish the antiferromagnetic ordering of energy 
levels in symmetric and antisymmetric sectors separately for wide range of the interaction constants. The nondegen- 
eracy of most levels is proved also. In Sec. Illli using the ordering rule and the results of Lieb and SchupjjiSiii on 
reflection-symmetric spin systems, the total spin and reflection quantum numbers of the ground state are derived. 
The results are tested and compared with numerical simulations. In Sec. lIVl the validity of obtained results is checked 
in some particular cases, for which exact results are already known. Section |V] is devoted to the frustrated ladder 
with boundary impurity spin. Similar ordering of energy levels is established in this case, too. Some partial results 
are obtained for the model with periodic boundary conditions in Sec. IVIl In the last section, we briefly summarize 
the results obtained in this paper. In Appendix, we apply the general results for reflection-symmetric spin systems 
obtained in Refs. [13 and [ll| to the frustrated ladder model. 



II. ENERGY LEVEL ORDERING OF FINITE FRUSTRATED LADDER WITH OPEN BOUNDARIES 



A. Frustrated ladder model 



In this paper we consider two finite identical antiferromagnetic spin-1/2 chains coupled by rung and diagonal 
interactionsi22i2i The Hamiltonian reads (see Fig.[T|): 




FIG. 1: Frustrated ladder with site-dependent couplings. Here jf are intrachain couplings while j/" and J;^ are correspondingly 
rung and diagonal interchain couplings. The dashed line is the symmetry axis of the model. 



N-1 N-1 N 

H = ^ jj (Si^i • Sij+i + 82,; • S2,;+i) + ^ J;^ (Si,i • 82,;+! + Sij+i • 82,;) + ^ Jt^i,i ■ ^2,;, (1) 
1=1 1=1 i=i 

where Si^; and 82,; are the spin operators of the first and second chains respectively. and are rung and 

diagonal interchain couplings, while jj' are intrachain couplings. The open boundary conditions are applied here. All 
couplings depend on site. Their choice corresponds to a model possessing the reflection symmetry TZ with respect to 
the longitudinal symmetry axis. The reflection just permutes the spins of two chains: 7?. 81,; — 82,(7?.. 

We do not put any restriction on the rung couplings Jj^ and consider antifcrromagnetic intrachain couplings only. 
We suppose also that the diagonal spin interactions are weaker than intrachain interactions, i.e., 

j}>\jr\- (2) 

Note that the condition above is not too restrictive for antiferromagnetic values of because the system remains 
invariant under the exchange of rung and diagonal couplings jf belonging to ^th box. The resulting model 

corresponds to a ladder obtained by permutation Si^p <-> 82,;' of two spins in all rungs positioned on the right from 
ith rung (i.e., for all V > I). It is topologically equivalent to the original model. Therefore, one can consider the 
couplings subjected to jf > only without loss of generality. Note that the condition ^ excludes the values of 

diagonal couplings equal to ijj'. 

One must mention that we call the model ^ a frustrated ladder, but, in fact, it becomes bipartite for some 
couplings. In particular, it is bipartite for ferromagnetic diagonal or rung interactions. These cases will be discussed 
in Sec. |TVl 

The Hamiltonian preserves the total spin 8 = Si^; + 82./ of the system. Since TZ and 8 are compatible, the 
eigenstates of ([T]) can be chosen to be parametrized by the spin {S = 0,1, . . . , N), spin projection (M = —N, —N + 
1, . . . , N), and reflection (cr = ±1) quantum numbers. 

In order to make easier the use of the reflection symmetry, we introduce the symmetric and antisymmetric super- 
positions of two spins on each rung: 

8|^) = 8i,/ + 82,,, S|°^ = Si,/-S2,/, l = l,2,...,N. (3) 

(s) . . . . . (s) 

The symmetrized spin 8j describes the total spin of ^th rung. It remains unchanged under reflection, i.e. 7?.8j TZ = 

sl'*^ while the antisymmetrized rung spin acquires minus sign, i.e., TZS'^'^^TZ = — 8[°\ 

Now we express the Hamiltonian ([T]) in terms of these operators. After omitting nonessential scalar term, the 
Hamiltonian takes the following simple form: 



N-l N 



Here, we have introduced the symmetrized and antisymmetrized couplings, which arc antifcrromagnetic due to the 
condition Q imposed above on the intrachain and diagonal couplings: 

tII , jx tII _ tX 

Jf ^ ' ; ' > and J- = ±—!^>0. (5) 
2 2 

Note that (j4]) does not contain terms, which mix symmetrized and antisymmetrized spin operators. This fact is a 
consequence of the reflection symmetry of the model. We mention also that similar decomposition for frustrated 
ladder Hamiltonian was applied also in Refs.HO and[2ll. 

In terms of lowering and rising operators S^''^'^ = S^^^^ ± iS^'^^^ and 5;"^^ = S^"'^'^ ± iS[°'''^ the Hamiltonian reads 



rj- _ 1 \^ ( js q{s)+ n(s)- . TSQ(s)-q{s)+ Taq(a)+q{a)- jaq{a)-Q 

n — - 2_^(Jl O;,^-^ -t- Jl + Jl + J; 



2 
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l+l > 



N-1 N 

1=1 



(6) 



Further, in this section, we will show that for jj > \J^\ the minimal energy levels Es.a in the symmetric (cr = 1) 
as well as in the antisymmetric (cr = —1) spin-5' sectors are nondegenerate and ordered antiferromagnetically, i.e.. 
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are increasing functions of S. This is an extension of Lieb-Mattis ordering theorem^ to the frustrated ladder model. 
Here, we briefly outline the steps of the proof. 

First, we construct a basis, in which all nonzero off-diagonal matrix elements of the Hamiltonian become negative. 
In the next step, we show that the matrix of Hamiltonian being restricted to the subspace of states with fixed values 
of spin projection = M and reflection TZ — a quantum numbers is connected. Then according to Perron-Frobenius 
theorem, the lowest-energy state in every such subspace, called a relative ground state, is nondegenerate. As we will 
show, the relative ground state in most cases has total spin value equal to absolute value of its z projection (S = \M\). 
Together with the total spin symmetry of the Hamiltonian this implies that the multiplet, to which the relative ground 
state belongs, possesses the lowest energy value among all spin-S* multiplets. The antiferromagnetic ordering between 
lowest energy levels Eg o- follows then directly from their nondegeneracy. 



B. Nonpositive basis 



The existence of a basis, where all off-diagonal elements of the Hamiltonian are nonpositive, is not obvious as it is 
for nonfrustrated models. The spin flip applied on one sublattic^»2§ is adopted for bipartite systems and does not lead 
to the desired basis for frustrated systems. In this subsection, we construct nonpositive basis for the frustrated ladder 
model. The first step is the use of the basis consisting of combined spin states on each rung instead of the standard 
Ising basis consisting of on-site spin-up and spin-down states. Then we apply a unitary shift to the Hamiltonian, 
which makes all nonvanishing off-diagonal elements arising from the symmetric part of ^ negative. In the final step, 
all basic states are multiplied by a sign factor making all nonvanishing off-diagonal matrix elements negative. 

It is rather difhcult to trace out the sign of matrix elements of H in the Ising basis, consisting of on-site spin states. 
Instead, we start with a more suitable basis consisting of combined spin states on every rung. Any rung state can be 
expressed as a superposition of three symmetric triplet states. 



|1>:=|1,1> = |J), |0):=|1,0) 







> 





and one antisymmetric singlet state, 



|0> 



(7) 



(8) 



We have used above standard notations "f, j for on-site spin-up and spin-down states. For convenience, we mark the 
triplet states shortly by labels ±1 and 0. 

The total space of states is spanned by the basic states 



mi) 



ImN) 



mi 



±1,0,0. 



(9) 



Here |m) is one of the four rung states defined above. The reflection operator TZ is diagonal in this basis. Its quantum 
number is (—1)-'^'', where Nq is the number of singlets in 

Now define unitary operator, which rotates the odd-rung spins around the z axis on angle tt, as follows: 



[(JV+l)/2] 



U = exp 



{s)z 
21-1 



(10) 



1=1 



where by [x] we have denoted the integer part of x. Under the action of f/, the odd-rung lowering-rising operators 

change the sign {Us!fi'_^l^U-^ = -Slfi'L'l^), while the others remain unchanged. Recall that for bipartite models, 
a similar unitary shift applied to the spins of one sublattice makes all off-diagonal elements of the Hamiltonian 
nonpositive In our case, the Hamiltonian © transforms to 



N-l 



1=1 

N~l N 



(11) 



1=1 



1=1 



It is easy to see that J"*" part of the Hamiltonian, presented by the last sum in (fTTj) . is diagonal in the basis ([9|) 
because it consists of the squares of rung spin operators. 
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The next observation is that J*^**) part of the Hamiltonian H gives rise to negative off-diagonal elements. Indeed, 
the symmetrized spin operators describe the spin of rung state and act separately on singlet and triplet states in the 
usual way. All their off-diagonal matrix elements are positive: 

(0| 1-1) = (1| |0) = |0) = (0| |1) = V2. (12) 

Taking into account the fact that the coefhcients J* are antiferromagnetic [see ([S])] , it is easy to see that all nonvanishing 
off-diagonal elements of the shifted Hamiltonian (jlip. which are generated by terms containing symmetrized spin 
operators, are negative. 

Finally, we consider the matrix elements produced by the antisymmetric local terms of Hamiltonian (|lip . In 
contrary to symmetric case, the antisymmetrized spin operators mix triplet and singlet states. All their nonzero 
matrix elements are off-diagonal and have the following values)^ 

(0|S'("'+ 1-1) = (-lis*'")- |0) = V2, (1|S'('^)+|0) = (0|5('^'-|1) =-\/2, (0| S-^")^ |0) = (0| S-^")^ |0) = 1. (13) 

Using the equations above, we can obtain all nontrivial matrix elements generated by J^"-* part of the Hamiltonian. 
The nontrivial action of terms with operators S*'"^)^ on two adjacent spins is 

^(a)T^(a)± 1^^^ ^ 1^^^ ^ _2 |0) ® |0) , S^^^^Si"-^^ |0) ® |0) = -2 |±1) ® iTl) , (14a) 

^(a)±^(a)T ^ 1^^^ ^ 2 |±1) ® |0) , ^^"^^5^"^^ |±1) |0) = 2 |0) ® |±1) . (14b) 

The action of terms in (fTT|) containing S*^"^^ reads 

s[^>si^^'\Q)®\Q)^\Q)®\Q), 5^"^"^^"^" |0) ® |6) = |0) ® |0) , (15a) 

5j"^"5^"^'|0)®|0) = |0)(»|0), 5^"^"5^"^"|0)®|0) = |0)®|0). (15b) 

The subscripts 1 and 2 indicate the rung, on which the spin operator acts. 

Taking into account the positivity of antisymmetrized couplings [see ([5])], we conclude that the elements 
generated by p4p acquire an overall minus sign in the shifted Hamiltonian matrix (jlip . At the same time, the 
elements corresponding to (fT5|) enter in the Hamiltonian with the same sign. Hence, only the contribution from (jl4bp 
gives rise to negative matrix elements in the Hamiltonian iJ, whereas (jl4ap and (|15p are responsible for unwanted 
positive off-diagonal elements. In order to alter their sign, we make the following observation. 

Due to the reflection symmetry, the parity of singlet states, i.e., (—1)^°, is conserved quantity under the action of 
each term in (fTTj) . The matrix elements ()14ap and (|15ap are the only ones, which are responsible for the creation and 
annihilation of a singlet pair. Other elements rest singlet number unchanged. Thus, multiplying the basic states ^ 
on sign factor (_i)™mber of singlet pairs ^ (_x)[Wo/2]^ g^e Can make the elements arising from (IT4a| and p5a| negative. 
The sign factor does not affect on other matrix elements. 

The remaining action (|15bp just permutes the singlet and 5^ = triplet states. In order to make this term 
negative, introduce the ordering between pair of rung states |0) and |0) inside multirung state dH). We say that a 
pair is ordered if |0) is located on the left-hand side from |0). Denote by N^^ the number of disordered pairs in ([5]), 
i.e., the number of pairs (|0) , |0)), where |0) is on the left-hand side from |0). Note that the nearest-neighbor actions 
(I15bp change the ordering of only one pair. Therefore, if we multiply a basic state by another sign factor (—1)^05, 
all matrix elements generated by (|15bp will change the sign and become negative. At the same time, other matrix 
elements will hold unchained. Indeed, it is easy to see, that the nearest-neighbor permutations |0) |±1) ^ |±1) |0) 
and |0) |±1) |±1) |0) do not change the number of disordered pairs TVqq while pair creations and annihilations 
(|±1) |=f1) |0) |0) and |±1) |=f1) |0) |0)) change its value on even number. Note that similar type of sign factor 
has been used in order to prove uniqueness of relative ground states for Heisenberg chains with higher symmetries 

Finally, the basis, in which all off-diagonal matrix elements of the Hamiltonian (llip are nonpositive, is 

|mi, m2, . . . ,mjv) := (-l)[^°/2]+JVog |^^^ ^ |^^^ . . . ® |mjv) • (16) 



C. Relative ground states in 5*^ = M, TZ — a subspaces 



Due to the spin projection and reflection symmetries, the Hamiltonian is invariant on each subspace with the definite 
values of spin projection and reflection operators: — M, TZ = cr, where M = —N, —N + 1, . . . ,N and a = ±1. 
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Keeping the terminology, we call it (M, a) subspace. Below, we outline the proof that the matrix of the Hamiltonian 
in basis being restricted on every (M, cr) subspace is connected. 

Note that due to ([5]) all local actions considered in (fT4|) . (jisp contribute in the Hamiltonian ([6|) with nonvanishing 
coefficients. It is easy to verify using ([T^ . p^ . and (fTS]) that any two adjacent states |m;) (g) |m;+i) and |toJ) (g) |toJ_|_j^) 
are connected by the ^th local terms of the Hamiltonian (jll[) provided that their quantum numbers are subjected 
to the conservation laws. In other words, both states must possess the same spin projection and reflection quantum 
numbers. This rule can be generalized by induction to any two basic states In fact, any symmetric (cr — 1) basic 
state (Uni) after successive applications of the local terms in pT|) can be transformed to state |±1,...,±1,6,...,0), 



where plus (minus) sign holds for positive (negative) values of M. Similarly, any antisymmetric (cr = —1) state is 
connected to state | ±1, . . . , ±1 , 0, . . . , 0, 0). This finishes the proof of the connectivity. 

|M| 

Now all conditions of Perron-Frobcnius theorem"^^ are fulfilled and one comes to the following result. 

• The relative ground state \^)Ma ^ *^ (M,a) subspace is unique and is a positive superposition of all basic 
states: 

> 0. (17) 

(_l)"o=<^ 



The state \^)m a- must have a definite value Sm,^ of total spin quantum number. Otherwise, it could be presented 
as a superposition of independent states with different spins but the same energy. This would be in contradiction with 
the uniqueness condition established above. It is evident that if some spin-5 state overlaps with the relative ground 
state \^) M a then S]\i^^ — S. Below, we use this property in order to determine Sm,it- 

Due to the spin reflection symmetry Sm.o- ~ S^M.a- So, one can consider only non-negative values of M. First, we 
suppose that M > 0. 

If CT = (—1)^"*^, then the state |0) — | 1, . . . , 1, 0, . . . , 0), which is the highest weight state for a spin S — M 

M N-M 

multiplet, contributes in the sum (fT7|) . Therefore, it overlaps with the relative ground state. According to the 
arguments above, Sm.o- — M for this case. 

Else if cr = (-l)^-^^-i, then both states |0, 1, . . . , 1, 0, . . . , 0) and |1, 0, 1, . . . , 1, 0, . . . , 0) are also presented in 

M N-M-1 M-1 N-M-1 

the decomposition of \^) jyj ^ with positive coefficients. According to the definition (jl6p of basic states, their sum, 

up to a nonessential sign factor, can be presented as {ip) = \/l/2 (|1) (g) |6) + |0) (g |1)) (g |1, . . . , 1, 0, . . . , 0) and, 
of course, overlaps with ()17p . In fact, has a definite spin value, which we will determine now. Remember 
we are currently working with the Hamiltonian (|lip . obtained from the original one by unitary shift (jlOp . Turn- 
ing back to H representation, we have to shift the states back. Thus, the original state is U^^ {ip) = U ~ 
\/T/2 (— |1) g) |0) + |0) g) |1)) g) 1 1, . . . , 1, 0, . . . , 0) up to an overall sign factor. The first term of the product is the 

M-l 

highest state of the triplet in the decomposition of two triplets, while the second term is the highest state of S* = M — 1 
multiplet. Together they form an S* = M multiplet. Therefore, Sm.o = M in the case of cr = (—1)^^^^^^, too. 

Consider now the two subspaces with M = 0. If cr = (—1)^, the state constructed from rung singlets participates in 
the sum pT|) . so the relative ground state is a spin singlet too: 5o,cr =0. If cr = (—1)^"^, the triplet state |0, 0, . . . , 0) 

enters in the sum, and, therefore, the relative ground state is a triplet, i.e., 5o,(j — 1. 
Summarizing the results, we arrive at the following conclusion: 

• The relative ground state \^)Mcr (always belongs to a spin-\M\ multiplet except M = and a = (—1)"'^^"'^ case 
when it belongs to a spin-triplet. In other words, 

5m. = r ^fM^Oanda^(-l)^-^ ^^^^ 
I ji ri other values of M and a. 
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D. Ordering rule among the lowest levels in sectors with different value of total spin 

Consider now the restriction of the Hamiltonian to the sector containing states with fixed values of the total spin 
and reflection quantum numbers. Apparently, the relative ground state \^) ]\j ^ is also the lowest-energy state in the 
sector characterized by 5* = ^m.ct and TZ = a. Denote by Es^a the lowest energy level in that sector. According to the 
previous subsection, the relative ground state is unique on (M, a) subspace. Therefore, the level Es,a with S = SM,<y 
must be nondegenerate, i.e., it must contain only one spin-S multiplet with parity a. This is just the multiplet, which 
encloses the state \^) m a itself. Using (fT8| . it is easy to see that any level Es.cr is nondegenerate in this regard, except 
perhaps the one with 5 = and a ~ (—1)^^^. 

Moreover, every (M, cr) subspace contains a representative from any multiplet with parity a and spin S > \M\. 
Therefore, the state \^) j^j ^ has the minimum energy for all these multiplets. According to p^ . its spin has the value 
SM,a = \M\ provided that a = (-1)^. The uniqueness of relative ground state then implies that all levels Es.a with 
spin S > Sm,(7 are higher than the level with S — Si\i,a-- Consequently, Es^a is a monotone increasing function of S. 
If cr = (—1)^"^, then the last equation is true for \M\ > only and, hence, Es.a increases in the range S > 1. 

Finally, we arrive at the following conclusion. 

• The minimum energy levels Es,a are nondegenerate (except perhaps the one with S ~ and a = (— and 



are ordered according to the rule: 



for a = ifSi>S2 

for a = (-1)^"^ if either Si> S2>1 or Si^ 0, 5*2 = 1. 



The ordering rule above enables to determine the total spin of the minimum-energy states in the symmetric and 
antisymmetric sectors. 

• The ground state in entire a ~ (—1)"'^ sector is a spin singlet while in a = (—1)^""'^ sector is a spin triplet. In 
both cases, it is unigue. 

In other words, for frustrated ladder with even number of rungs, the ground state in the symmetric sector is a 
singlet state, while in the antisymmetric sector, it is formed by the three states of a triplet. In contrary, for odd 
number of rungs, it is a singlet in the antisymmetric sector, while in the symmetric sector, it consists of three triplet 
states. 



III. GROUND STATE 



A. Exact results 



The ordering rule (|19p does not compare energy levels in symmetric (cr = 1) sector with the levels in antisymmetric 
(cr — —1) sector. It is clear that the total ground state coincides with the minimum-energy state in either symmetric 
or antisymmetric sector. If the lowest levels in both sectors coincide, the total ground state becomes degenerate. 
Using also the results of previous section, we come to the following conclusion. 

• The ground state of frustrated ladder with couplings obeying ^ may be: (a) a unique a = (—1)^ spin singlet; 
(b) a unique a = (— 1)^"-'^ spin triplet; (c) any superposition of both of them, i.e., singlet® triplet. 

The frustrated ladder ([I]) in certain parameter region belongs to a large class of reflection-symmetric models consid- 
ered recently by Lieb and Schupp.^^ In Appendix, the rigorous results obtained in Ref . 1 1 1| are applied to the frustrated 
spin ladder case. It appears that all ground states of the Hamiltonian ([T]) are spin singlets if the couplings satisfy 
Jj^ > \Ji^i \ + \Ji^\, where Jq = = is supposed. The nonstrict inequality sign can be set instead of strict one 
for antiferromagnetic values of diagonal couplings. The details are given in Appendix. Note that the aforementioned 

result is true for any values of intrachain couplings jf including ferromagnetic ones. 

From the other side, for jf > |J;^| our results established in previous section become valid too. Setting together 
both Lieb-Schupp and our results, we come to the conclusion that only case (a) above can take place. This means 
that the ground state is a unique cr = (—1)^ singlet, provided that the couplings satisfy 

Ji^>\jr-i\ + \jn and jf>\jn- 

Consider separately the most familiar case of ladder with site-independent couplings. According to the discussions 
above and in Appendix, we come to the following conclusion. 
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FIG. 2: Two lowest energy levels, obtained by exact diag- 
onalization of frustrated ladder with 18 sites (A'^ = 9), are 
plotted as a function of rung coupling. Corresponding states 
have different spin S and reflection a quantum numbers. The 
level crossing happens at some point. 



FIG. 3: The ground state and lowest excitation of frustrated 
ladder with 20 sites (A'^ = 10) obtained by exact diagonaliza- 
tion for different values of rung coupling. Both states have the 
same quantum numbers. The level crossing does not happen 
in agreement with our exact results. 



• The ground state of finite size frustrated ladder model with antiferromagnetic rung and intrachain inter- 
actions is a unique spin singlet with a = (—1)^ reflection symmetry if the couplings satisfy > 2J^ > — J^ 
and Jll > |J^|. 

It must be emphasized that this result is rigorous for finite size ladders only, which are the main objects of study 
in this paper. 



B. Thermodynamic limit, comparison with other approaches and numerical test 

The results obtained above are in a good agreement with the ground-state properties of frustrated ladder model, 
which have been investigated intensively in the literature so far. From our results, it is unclear whether or not the 
singlet-triplet degeneracy takes place for finite-size ladders. Moreover, in the thermodynamic limit oo, the 

additional degeneracy of energy levels can occur and, in principle, the strict inequality in p9p must be replaced by 
nonstrict one. The investigations carried out by different methods suggest that all three possibilities for the ground 
state described in the previous subsection may take place. In fact, the ground-state degeneracy happens at critical 
points only, while in most cases, the ground state is a unique a = (—1)^ singlet or cr = (—1)^"^ triplet. This property 
remains true in the thermodynamic limit. The rest of this section is devoted to the testing of our exact results and 
more detailed comparison with data obtained by other methods. 

The weak interchain coupling analysis (J^, ^ j") based on conformal field theory approach shows that for 
> 2J^, the ground state has a tendency to form singlets along the rungs and triplets along the diagonals, while 
for < 2 , it has a tendency to form triplets along the rungs and singlets along the diagonals. The former ground 
state corresponds to the usual ladder or rung-dimer phase, while the last one corresponds to spin-1 Haldane chain. 
In both cases, it is nondegenerate. The two phases belong to the same universality class because a continuous path 
connecting them exists and it does not contain any critical point. "^^ The distinction between them has a topological 
character ^?^''^M For = 2J^, both types of dimerized ground states exist leading to twofold degeneracy. ^'^ 

Numerical simulations2ii2& strongly support this picture and extend the phase transition curve out of the weak 
coupling region. The transition curve obeys the relation < 2J^ and approaches to the line = 2J^ at the weak 
coupling limit. 

The properties of the ground state differ for ladders with odd and even values of the rung number. First, we consider 
odd-rung ladders. Small values of rung coupling give rise to Haldane-type ground state, which is a spin triplet 
with a — 1 [case (b) in the previous subsection] while large values lead to rung-singlet-type ground state with a — — I 
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[case (a)]. At some intermediate point, the level crossing happens, which corresponds to aforementioned case (c) with 
degenerate singlet-triplet ground state. The level crossing point approaches to the critical line in the thermodynamic 
limit. The numerical investigation of small size ladder Hamiltonians supports this picture. In Fig. [21 the two lowest 
levels of frustrated ladder with 18 sites {N — 9) obtained by exact diagonalization are plotted for different values of 
rung coupling. The figure clearly indicates that both levels cross at some point. 

On the contrary, for even-rung ladders, both Haldane and rung-dimmer states are singlets and located in cr = 1 
sector. The level crossing is forbidden in this case because our results exclude the existence of degenerate double-singlet 
ground state for finite-size systems. In Fig. [31 we present the two lowest levels obtained by numerical diagonalization 
of the system with 20 sites {N = 10) close to the composite spin point j" = . Both levels approach each other but 
do not intersect. The level crossing is absent. It occurs at the composite spin point, which is out of the parameter 
space ^ considered in this paper. This model will be discussed in the next section. In the thermodynamic limit, the 
closest point of both levels approaches to the level crossing point of odd-rung ladders and the level crossing picture 
recovers. 

Similar behavior is observed for ferromagnetic interchain couplings (J^ < and < 0). Large values of \J^\ 
correspond to Haldane phase and its small values put the system in rung-dimmer phase. The phase transition curve 
in the thermodynamic limit is remarkably close to the line = 2J^ (Ref. 33) and coincides with it in the weak 
interchain interaction limit.— This curve has been traced out with high precision by applying different boundary 
conditions for each phasei^ A slight deviation from it is observed close to the ferromagnetic phase only. 



IV. COMPARISON WITH KNOWN RESULTS 



In this section, we compare the properties of the generalized ladder Hamiltonian (|T]), ^ established in the previous 
section with those obtained earlier for some particular values of couplings. 



A. Ferromagnetic rung interactions 

If we set all rung couplings to be ferromagnetic and others to be antiferromagnetic, the ladder model ([T]) becomes 
bipartite. The two sublattices A and B, forming the ladder, consist of spins of even and odd rungs correspondingly. 
Then the conditions of standard Lieb-Mattis ordering theorem^ are fulfilled. The Ising basis after the unitary shift 
(ITUl) becomes negative. Thus the relative ground state is unique in whole subspace = M. Apparently, 

it is an eigenstate of the reflection operator. The corresponding eigenvalue is cr = 1 because |f2)^,j is a positive 
superposition of shifted basic states^ and the shift operator (|10p remains invariant with respect to the reflection. 
Thus, in this case the minimal energy in the symmetric sector is lower than the corresponding one in antisymmetric 
sector: Es — Es^i < Es^-i- As was mentioned in the Introduction, the ground state is unique and has spin 
Sgs — \Sa — Sslr^ i.e., is a singlet for even values of N and triplet for odd values of N. In both cases, a = 1. This 
is in agreement with our results, which claim that a = (—1)^ for singlet ground state and a = (—1)^^-'^ for triplet 
one. All mentioned properties are similar to the properties of spin-1 Haldane chain because due to ferromagnetic 
rung couplings, which make preferable rung triplet states, the low-energy properties are well described by Haldane 
chaini^i^^ 



B. Ferromagnetic diagonal interactions 



Consider now the ladder with ferromagnetic diagonal couplings and antiferromagnetic rung and intrachain couplings. 
The low-energy properties of this system have been investigated recently in Ref. 33. The model contains standard 
antiferromagnetic ladder as a particular case. The frustration is lost in this case, too. The lattice becomes bipartite 
of checkerboard type and the Lieb-Mattis theorem can be applied also. The sublattice A consists of the odd sites 
of the first chain and the even sites of the second chain, while the sublattice B consists of the even sites of the first 
chain and the odd sites of the second one. The spin-rotation operator leading to a nonpositive basis is given by 
U' = exp {iirS^), where S\ = X^zC-^i 2i-i + ^221) ^Pi"^ projection operator of A sublatticci^ As in the previous 

case, the relative ground state |r2)^ in M subspace is expressed as a positive superposition of the shifted basic states 

where m =1, | labels on-site spins. Now using the definition of TZ and U' 
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FIG. 4: Spin-1/2 frustrated ladder with boundary impurity spin sq. 

Hence, we have TZ \^) j^j = (—1)^^^'^ Thus, among two (M, a) subspaces with different values of cr, the minimal 

energy level in the subspace with cr = (—1)^^*^ is the lowest: Es — Es.a < Es-a- The total ground state is a spin 
singlet with a — (—1)^ reflection quantum munber in accordance with our results in Sec. |lll 

C. Composite spin model 

The equal values of diagonal and intrachain couplings {J^ = jf) correspond to so-called composite spin modeli^ 
These values of coupling parameters are out of the range ([2]), where the results of Sec. |n]are valid. However, due to 
continuity, some results still remain true in this limiting case. In particular, the ordering rule is fulfilled after replacing 
the strict inequality sign in (|19p by the nonstrict one like in the thermodynamic limit. 

In contrast, the results related to the uniqueness of relative ground states cannot be applied here. Indeed, in this 
case, Ji'^^ — and the antisymmetrized term, which is responsible for the exchanges between singlet and triplet rung 
states, is absent in the Hamiltonian (llip . The spins of individual rungs are now conserved quantities.-^" Therefore, 
all singlets remain frozen at their sites and the Hamiltonian is not connected on (M, a) subspaces anymore. For 
site-independent values of couplings, the ground state has been obtained explicitly^ For small values of J-"-, it is 
reduced to the ground state of Haldane chain (5 = for even N, S = 1 for odd N and cr = 1 always), while the large 
values lead to the rung-sing let state |0) (g) • • • (g) |0) [5 = and cr = (-1)^]. This agrees with our results in Sec. IIIII 
Both states are eigenstates of the composite spin Hamiltonian, which makes the level crossing at some critical point 
J^{N) < 2 inevitable. For even values of N, the degeneracy disappears under a weak deviation from the composite 
spin point, as has been discussed in Sec. IIIII However, the two lowest levels still approach closely, as is shown in 
Fig. [31 

V. ORDERING IN FRUSTRATED LADDER WITH BOUNDARY IMPURITY 

Consider frustrated spin-1/2 ladder ([T]) with an impurity spin sq coupled antiferromagnctically to two spins of the 

first rung, as was shown on Fig. [3) The original Hamiltonian ([T]) is supplemented by boundary term H\aip = Jo Sq-S^*', 
Jo > 0, which preserves the reflection symmetry. 

A. Relative ground states 

After the rotation of odd-rung spins on angle tt around z axis by means of unitary operator (|10p , the impurity-ladder 
interaction term acquires the form. 

Hirap = UHi^pU = ^{^o'^l ^ -\- Sq S^^ ) + JqSqS^ ^ . 

We will work here in a basis, which is a natural extension of the basis (ITB|) . used before for the ladder with free 
boundary conditions: 

|too) (g |mi, . . . , itin) = (-1) [^0/21+^00 \mo) |mi) (g . . . (g Itun) ■ (20) 

Here |mo) := |so,m.o) (toq = —sq,...,sq — 1,sq) is the usual basis of spin-so multiplet. All nonvanishing matrix 
elements of Sq are positive: (mo| Sq \m.Q — 1) = (too — 1| |too) = (so + mo) (so — Too + 1) > 0. It is easy to see 
that the spin exchanges due to the boundary term do not affect on the sign factor (-l)[^o/2]+^oo. This together with 
(fT^ implies that Jfimp has only nonpositive off-diagonal elements in the basis (PO)) . 
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The bulk part (jlip of the Hamiltonian is also nonpositive in this basis as was proved already in Sec. |TT1 Thus, the 
entire Hamiltonian H + Himp has no positive ofF-diagonal elements, too. 

It is easy to ensure that its restriction to any (M, a) subspace with fixed spin projection M and reflection a quantum 
numbers is connected. This can be shown using the arguments similar to those applied before for the ladder without 
impurity spin. 

Applying again Perron-Frobenius theorem, we come to the conclusion that the relative ground state of i7 + Hin\p 
of each (Af, a) subspace is unique and is a positive superposition of all basic states (|20|) : 

>0. (21) 

(-1)"0=^ 

Now turn to the determination of the spin S'j\/,cr of this state. As before, in order to do this, we will construct a 
state with certain value of spin and positive or vanishing coefficients in its decomposition over the basic states ()20|) . 
In this case, it will overlap with the relative ground state (PT|) . Then the uniqueness of the relative ground state would 
indicate that both states have the same spin value. 

Due to the spin reflection symmetry, both states have the same spin quantum number. Thus, without 

any restriction, we can consider M > values only. First, we consider the values M > sq and will prove below that 

SM,a = M. 

In a ~ (^_i^N-M+so sector, there is a simple highest state |so) (8) | 1, . . . , 1, 0, . . . , 0) of a multiplet with spin S ~ M. 

M-so 

It coincides up to an unessential sign factor with one of the basic states d^U]) . Therefore, Smm — M for this case. 
For the opposite values of the reflection quantum mmrber [<t = (— l)^~*^+*o-i]^ -^g ^g^j^g g^g a,n indicator the state 

1^') = |.so)® |6) + ^=i== |.o-i>®|i>) ®|i^_^,o,...,o), 

A/-SO 

which can be considered as a generalization of the state \ip) used before for bulk ladder. Again, \ip') is a positive 
superposition of two basic states from ([20|) . As has been done before for state lip), its spin can be detected going 
back to the original H representation: lip') |'0'). As a result, just the sign of the second term in the brackets 

is changed. It is easy to check that \so) ® |6) — ko — 1) <8i |1) is the highest state of spin-so multiplet, 

which appears in the decomposition of the tensor product of two multiplcts with spins sq and 1. Thus, the spin of 
the entire shifted state \ip') is again M. 

Next, we consider the values < M < sq, which are essential for sq > 1- We will show that in this case the relative 
ground state in {M,a) subspace belongs to a spin SM,a = sq multiplet if cr = (—1)^; otherwise it belongs to a spin 
SM,a = So — 1 multiplet. The first case is easy to prove using for the test the basic state \M) (g) |0, . . . , 0) of spin sq. 

For the second case we choose as a test the state \ip") = |x) ^ |0) ® . . . ® |0), where 

^ ^ ^ 



'(.so- 


-A/+l)(.so- 


M) 


2so(2so + l) 




f M+l)(so4 


-M) 




|M + l)(g) 1-1). 



2so(2so + l) 

It is easy to see that IV'") is a positive superposition of three states from the basic set ([20| . The action of on 
IV'") changes only the sign of the first and last terms of \x) giving rise to a state of spin 5 = Sq — 1 multiplet, as is 
easy to verify by applying (So + S^''^)^ to it. 

Finally, summarizing the results obtained above, we arrive at the following conclusion. 

• For the ladder with diagonal interactions obeying ^ and antiferromagnetically coupled boundary impurity sq 
the relative ground state in {M, a) subspace is nondegenerate and belongs to a multiplet with the spin 

*/ \M\ > so 

SM,a ^ {so if \M\ < So and a = (-1)^ (22) 

1 i/|Af| < So and cr = (-1)^-1. 



Note that for so = 1/2 value of impurity spin, only the first line holds and we have simply SM,a- = \M\. 
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B. Ordering rule and ground state 

Using the spin value of the relative ground state (|22p and its uniqueness proved above, one can compare the minimal 
energy levels Es^a- of the sectors with fixed spin and reflection quantum numbers and study their degeneracy just in 
the same way as was done above for the open ladder. 

For (T = (—1)^, it is clear from (|22p that the energy level Es.a is nondegenerate for any S > sq. The uniqueness 
of the relative ground state suggests that all levels with spin S > Sm,^ are above the level with S — SM,a, which 
contains the state \^) m a itself- Thus, Es^a is a monotone increasing function of spin in the range S > sq- For S < sq, 
the level degeneracy and ordering are not clear. All relative ground states \^)m a ^i^h \M\ < sq constitute a spin-so 
multiplet having the lowest-energy value among all states with the reflection parity a = (—1)^. 

For a = (—1)''^"^, the level Es^a is nondegenerate for any S > sq — 1. The lowest energy Es^a is a monotone 
increasing function of spin in the range S > Sq — 1. In the case of sq > 1, all relative ground states with 

\M\ < So — 1 are combined into S — so — I multiplet having the minimal energy among the states with a = (—1)^^-'^. 
For So = 1/2, the lowest level increases with spin everywhere. The following statement sums up the discussions above. 

• For the ladder with diagonal interactions obeying jf > | Jj^ | and antiferromagnetically coupled boundary impurity 
So the lowest energy levels in sectors with fixed value of spin S and reflection cr quantum numbers are ordered 
according to the rule: 

(-1)^ if either Si > S2 > So or Si < S2 ^ So ^^3) 
(— l)^^""^ if either Si > S2 > so — I or Si < S2 ~ so — 1- 

• The ground state of the model in entire a = (—1)^ sector is a S = so multiplet, while in entire a = (—1)^"^ 
sector is a S" = |so — 1| multiplet. In both cases, it is unique. 

The absolute value sign is used in order to take into consideration the impurity with sq — 1/2 also. In that case, 
the ordering rule ([23)) is simplified because in that case S — 1/2 is the lowest value of the total spin and the relation 
Esi.a > Es2,a- holds for any a and any two spins obeying Si > S2. The lowest-energy states in both symmetric and 
antisymmetric sectors are unique spin doublets. 

As for the ladder model without impurity, one cannot compare energy levels related to different reflection quantum 
numbers. Nevertheless, one can gain information about the total ground state of the model from the ordering rule 
established above. The total ground state is nondegenerate and cither belongs to cr = (—1)^ sector and, hence, has 
spin value 5 = so or belongs to cr = (—1)^^^ and has spin = |so — 1|. If the lowest energy levels in both sectors 
coincide, which can occur, for example, due to some additional symmetry presented in the model, then the ground 
state becomes degenerate and is the superposition of two multiplets with spin values sq and jsp — 1|. Therefore, the 
ground state could be at most doubly degenerate. 

In contrast to the model with free boundary conditions studied in Sec. |TT1 here we cannot apply the Lieb-Schupp 
resultsiSiiii in order to choose the valid quantum numbers for the ground state. The reason is that their approach 
cannot be used for a reflection-symmetric spin system if some spins are positioned on the symmetry axis. 



Esi,<j > Es2,a ■ I 



C. An example: Periodic spin-1/2 chain with odd number of spins 

The translationally invariant periodic spin-1/2 chain with odd N' = 2N + 1 number of spins is a particular case of 
the general class of frustrated ladders with boundary spin, considered in this section. It corresponds to sq = 1/2 and 
Jo = jf = Jjv with vanishing values for other coupling coefhcients. The model is integrable by Bethe ansata^''' like 
its more familiar even-site counterpart, but many of its properties differ from those of even-site chaini^ In particular, 
the even-site chain is bipartite and Lieb-Mattis theorem is valid in this case, while the odd-site chain is a frustrated 
system. Recently, the classical ground state of the last model has been constructed and the lowest energy has been 
obtained exactlyi^S 

The system possesses an additional symmetry: it remains invariant with respect to the cyclic translation by one 
site T. The translation and reflection operators satisfy the commutation relation TZT = T~^TZ. Thus, if is an 
eigenstate of T with eigenvalue e*"^ then the reflected state TZ is also an eigenstate with eigenvalue e~^'^. Hence, 
for all values of momentum except (/< = 0, tt, the energy levels of periodic chain are at least twofold degenerate. 

The relative ground states of even-site translationally invariant chain have just the exceptional values of momenta (0 
or tt), which is in agreement with their uniqueness. In contrast, for odd-site chain, the exact solution shows that relative 
ground state in 5*^ = M subspace is exactly doubly degenerate with two opposite momenta <j) — ±7r(M -I- l/N')M- 
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This is in agreement with our results, which assert that it should be at most doubly degenerate. Therefore, the 
relative ground states in symmetric and antisymmetric M subspaces have the same energy levels. The reason of 
the degeneracy is the translation symmetry, which mixes these two states. Hence, the lowest levels in symmetric 
and antisymmetric sectors coincide (i5s,i = Es.-i — Es) except, of course, M — ±7V'/2 subspaces, where only one 
ferromagnetic state presents. Moreover, antiferromagnetic ordering rule {Es^ > i?s2 if > S2) is fulfilled like for 
even-site chains. The total ground state is the combination of two spin doublets with different reflection quantum 
numbers. 



VI. PERIODIC LADDER 



Consider now the frustrated ladder with periodic boundary conditions. The bulk Hamiltonian ([T]) is supplemented 
by the term 

^^por = J'Jv(Si,Ar • Sij + S2,Ar • 82,1) + Jj^(Si_jv ' S24 + Sij • S2,Ar), 

which describes the interactions between the spins of the first and last rungs. Following ([2]) we set the condition 
> I I on the boundary coupling constants. In the expression of the Hamiltonian, in terms of symmetrized and 

antisymmetrized spin operators (g]), the boundary term contributes as Hp^r = Jn S^'' ■ sj*' + sj^^ • S^"\ where the 
coefficients J^, are given in Both are antiferromagnetic. 

Under the action of unitary shift operator (|10p . the total Hamiltonian acquires the form H + -ffpcr, where H is 
derived already in pT|) and 



(24) 



The boundary term iJpor produces positive off-diagonal elements in the basis (jl6p . in which the bulk part H of 
the total Hamiltonian is negative. One reason is the sign factor (—1)^^^ in ((24|) due to which positive off-diagonal 
elements appear for even values of N. Another more important reason is that the sign factor {—l)^oo in front of basic 
states p6l) is essentially nonlocal. Remember that JVqq counts the number of such (|0) , |0)) pairs that |0) is positioned 
on the left-hand side from |0). The exchange between the first and last rung spins can produce an uncertain sign 
factor, which depends on intermediate rung states. Below, we will derive the conditions, under which the basis (|16|) 
nevertheless remains nonpositive. 

The term S'j;''''^iS'j^-'^ is diagonal in the basis (|16p . The other terms produce off-diagonal elements. First, consider 
the matrix elements, which are generated by the terms with 5^°'^. They do not depend on the parity of N. Using 
(fT5)) and pB]) it is easy to check that 

(6, . . . , 0| s[^^'S^j^'>' |0, . . . , 0) = (0, . . . , 0| si^^'S^j^^' |0, . . . , 0) = (-1)^0+^0+1, (25) 

where the dots indicate all intermediate rung states, which have to be the same for the bra and ket states to ensure 
that the matrix element is nonzero. Here, Nq and Nq are the numbers of rung-singlet and — triplet states 
correspondingly. 

Next, we consider the matrix elements generated by lowering- rising spin operators in (I24p . Taking into account 
(IT2|) . (fT4l) . and we obtain 



(6, . . . , 6| si'^^s^^^^ |±i, . . . , Ti) = (±1, . . . , o| si'^^s'^^^^ |0, . . . , ±1) = (-1)^" • 2, 

(0, . . . , 0| si^'^^S^j^^^ |±1, . . . , Tl) = (±1, . . . , 0| si^'^^sP^ |0, . . . , ±1) = (-1)^6 • 2. 

The matrix elements (|25|) and (|26|l together with their transpositions are the only off-diagonal elements generated by 
iJpor. They will be negative if (-l)^o+^o = (_i)A'+JVo = (_i)JV+7Vo ^ i because the values of coupUngs and 
are positive. Using the relations <j — (—1)^" for reflection quantum number a and (—1)^^ = (— l)^~^fl~^o fQj- ^j^g 
spin projection quantum number M , we can rewrite the last equations as 

a=(-l)^ = (-l)^'^. (27) 

Remember now that the bulk part of the Hamiltonian has only nonpositive off-diagonal elements in the basis (fT6|) 
for any M and a (see Sec. [Hi. Therefore, the entire Hamiltonian is nonpositive only on those (M, cr) subspaces 
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(where = M and TZ = a), which are subjected to the condition ((27)) . Moreover, in any such subspace, its matrix 
is connected because H is connected as was aheady proven. Thus, according to Perron-Frobenius theorem, the 
relative ground state in each (M, a) subspace with quantum numbers obeying ([27]) is nondegenerate and is a positive 
superposition of all basic states (fT7|) . The last fact has been used in Sec. [IT] in order to determine the spin S'm.o- of the 
relative ground state for the open ladder. So, here, the formula (jlSp can be applied too, and, actually, it is reduced 
just to SM,a — \M\ since the equations (|27|) exclude the exceptional case of Af = and a — (—1)^"^. The outcome 
is as follows. 

• For the periodic frustrated ladder with even number of rungs, the relative ground state in any (M = even, cr = 1 j 
subspace is unique and belongs to a spin-\M\ multiplet. For the ladder with odd number of rungs, the same is 
true in any (M = odd, a = —1) subspace. 

As a consequence, the lowest energy levels i^s^o- among all states with fixed spin S and reflection a quantum numbers 
are nondegenerate for a = I and even values of S as well as for cr = — 1 and odd values of S. Moreover, we get a 
partial antiferromagnetic ordering of the lowest energy levels Eg o-- 

• For the periodic frustrated ladder with even number of rungs, the relation i?Si,i > Es2.i holds for any even spin 
S2 and any spin Si such that Si > S'2. The ground state in entire symmetric (a = 1) sector is a unique spin 
singlet. 

• For ladder with odd number of rungs, £^Si,-i > ^$2,-1 for any odd spin S2 and any spin Si satisfying Si > S2. 

For odd N, the lowest state in the antisymmetric may be a unique spin triplet, a spin singlet(s), or a 
singlet(s)© triplet representation. 

Remember that in this paper we consider the ladders with coupling obeying jf > \J{^\. The ladder with periodic 
boundary term under certain supplementary condition on coupling constants becomes a part of more general class of 
reflection-symmetric models investigated recently by Lieb and Schupp.^"'^^ This fact has been used in Sec. [IT] in the 
case of free boundaries and is discussed in detail in Appendix. For the coupling values obeying J^^ > \Ji^_i \ + | J;^ |, all 
ground states become spin singlets, and among them there is one with a = (—1)^. Moreover, for antiferromagnetic 
diagonal interactions (i.e., if all > 0) the nonstrict inequality sign can be used in the last inequality instead of strict 
one [see (jA7[) in Appendix]. In particular, the translationally invariant frustrated ladder with the antiferromagnetic 
couplings obeying J-^ > 2 and j" > has only spin-singlet ground states as well as is subjected to the ordering 
rule proved in this section. 

For the ladders with even rungs, the aforementioned properties are in good agreement with our results obtained in 
this section. For the ladders with odd rungs, they even strengthen the partial ordering rule. The relation Es^.-i > 
Es2,-i, holding for any odd S'2 and any Si such that Si > S2, becomes valid for S'2 = 0, too. The question of the 
degeneracy of the level -Bq.-i still remains open. 

Unlike in open ladder case, here we did not obtain definite exact results related to the degeneracy degree and spin 
quantum number of the total ground state. The main reason is the absence of exact results related to cr = (—1)^^^ 
sector in the case when the periodic boundary term presents. 



VII. SUMMARY AND CONCLUSION 



We have generalized Lieb-Mattis energy level ordering rule for spin systems on bipartite lattices^ to the frustrated 
spin-1/2 ladder model. The model consists of two coupled antiferromagnetic chains frustrated by diagonal interactions 
and possesses reflection symmetry with respect to the longitudinal axis. The spin exchange coupling constants depend 
on site and are subjected to the relations We have considered the finite-size system with free boundaries, the 
system with any impurity spin attached to one boundary as well as the ladder with periodic boundary conditions. 
Below, we describe briefly the results obtained in this paper. 

The total spin S and reflection parity a = ±1 are good quantum numbers. So, the Hamiltonian remains invariant 
on individual sectors with fixed values of both quantum numbers. For the open ladder and ladder with impurity, we 
have established that the lowest-energy levels Es.a of these sectors with the same value of a are ordered antiferro- 
magnetically. More precisely, the relation Es^^a > Es2,a holds for any two spins satisfying Si > S2 > Sgs(cr). Here, 
Sgs((T) is the spin value of the unique lowest level multiplet among all multiplets with the reflection parity a: 

. , , f^o ifa = (-l)^ 

"^^^''^-ll.so-ll if a = (-1)^-1, 
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where N is the number of rungs and sq is the impurity spin value. The case Sq = corresponds to the open ladder 
without impurity. We have also proven that all levels Eg^^ with S > Sgsi<^) are nondegenerate, which means that 
only one multiplet exists on the corresponding level. Note that S'gs (cr) coincides for (7 = 1 with the ground-state spin 
Sgs of the Haldain chain obtained after replacing each rung with the spin s = 1. The last model is bipartite and the 
Lieb-Mattis formula Sgs — \Sa — Sb\, mentioned in the Introduction, gives just the value of S'gs(l). 

In contrast to bipartite spin systems, the lowest levels in the symmetric (cr = 1) and antisymmetric {a = —1) sectors 
cannot be compared with each other, at least by means of our approach. Therefore, in general, the total ground state 
of the model is a unique multiplet with spin sq or |so — 1|. If the lowest levels in both sectors coincide, then the 
degeneracy occurs and the ground state is the superposition of both multiplets. In this regard, our results claim that 
the ground state can be at most doubly degenerate. 

For more restrictive values of couplings given in Appendix, the frustrated ladder with free boundaries (sq = 0) fits 
the class of reflection-symmetric spin systems, all ground states of which are spin singlets J^i^i Combining with our 
results, this implies that the ground state is a unique spin singlet in this case. This property is true, in particular, for 
frustrated ladder with site-independent antiferromagnetic couplings obeying J-^ >2J^ and j" > J^. 

For ladder with periodic boundary conditions, we have proven a weaker ordering rule. Namely, for the ladder with 
even (odd) number of rungs the ordering Es-^.a > Es2.a if Si > S2 is established exactly for a = 1 {a = —\) and even 
(odd) values of spin S2 only. The degeneracy of the ground state and its total spin value remain open questions in 
this case. 
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APPENDIX A: APPLICATION OF LIEB-SCHUPP APPROACH TO FRUSTRATED SPIN LADDER 

Recently, Lieb and Schuppii proved exactly that a reflection-symmetric spin system with antiferromagnetic crossing 
bonds has at least one spin-singlet ground state. Moreover, under certain additional conditions, all ground states 
become singlets. Frustrated spin-1/2 ladder ^ is an example of reflection-symmetric system (see Fig. [Ij. The 
Hamiltonian in the Lieb-Schupp form is 

AT JV N 

H li^i.i ■ ^2,1 + ^{aiSij + PiSi^i+i) ■ {aiS2,i + AS2,i+i) + ^ M^i.i ■ ^i.i+i + 82,; • S2,i+i). (Al) 
1=1 1=1 1=1 

Here, ai, f3i, 7;, Ji are arbitrary real coefficients. The boundary conditions are ckat = (3n = Jn = for open ladder and 
S5,Ar+i = S^^i ((5 = 1, 2) for periodic ladder. There are two kinds of bond in (|A1[) : the first sum consists of rung bonds 
while the second sum contains square bonds. Comparing (|A1[) with ([1]) one can express the spin coupling constants 
in ([1]) in terms of Lieb-Schupp coefficients: 

jll = Ji , Jt ^af+ Pi 1 + 7f , = A • ( A2) 

In the second equation above the boundary conditions /3o = and /3o = Pn are used for open and periodic ladders 
correspondingly. Note that the value of interchain rung coupling is always antiferromagnetic. 

Exploiting the reflection symmetry, Lieb and Schupp proved that there is a ground state of (|A1[) . which overlaps 
with rung-singlet state |0, . . . , 0) (Ref.[l3) and, more generally, with the canonical singlet statej^i This means that the 
spin ladder Hamiltonian with diagonal interactions ([T|) possesses a spin-singlet ground state with a = (—1)^ reflection 
symmetry for any value of jf provided that the remaining two couplings are subjected to the relation 

Ji^>af+Pli, J{=aiPi (A3) 

for some ai, Pi. Setting several restrictions on auxiliary parameters ai,Pi,ji, one can simplify (|A2[) and (jA3[) signif- 
icantly and get the constraints containing interaction couplings only. Of course, those restrictions narrow the total 
region in phase space, where Lieb-Schupp result is valid. 

Consider the simplest and most familiar case of ladder with site-independent couplings. We set ai = a, Pi = P, 
7i = for periodic ladder, where a, P are some real parameters. For open ladder, we choose the boundary values of 



16 



7/ parameter as 71 = P,"fN = ct. It is easy to verify that both for open and periodic boundary conditions, all values 
of interchain couphngs satisfying J-^ = + /3^, = a/J lie inside the general region in phase space defined by (|A3p . 
Hence, there is a ct = (—1)^ singlet among the ground states if two interchain coupling constants satisfy 

J^>2|J^|. (A4) 

This relation can be generalized to the ladder model with site-dependent couplings. Setting ai — ±/3;, one reduces 
the general region (jA3|) to 

Ji^>\jr-^\ + \jr\- (A5) 

Here, as before, the boundary conditions Jq = = and Jq ~ are applied for open and periodic ladders 
correspondingly. 

Another important property of reflection-symmetric spin systems established exactly in Refs. [l^ and fill is that for 
any crossing bond, the expectation of the spin of its all sites, weighted by their coefficients, vanishes for any ground 
state |ri). Using this feature of ground states, it was proved that in a system with sufficient symmetry, so that every 
spin can be considered to be involved in a crossing bond, all ground states are spin singlets. 

As was already mentioned, for frustrated ladder, there are rung bonds and square bonds, which are described 
correspondingly by the first and second term in (jAip . The aforementioned condition means that {n\ sl"^ |r2) if 
7j 7^ and ajSl"-* + PiS'l'^^ \ — 0, where sj"'' is the spin operator of Zth rung ([3]). Therefore, if a; = /3; 7^ 0, 

then the expectation of total spin of each square bond vanishes also, i.e., s[''-' + S^'^^ \n) — 0. Applying this to the 
ladder with site-independent spin exchange, one can ensure that all its ground states are singlets provided that the 
interchain couplings satisfy 

> 2J^ > -J^. (A6) 

Indeed, using aforementioned parametrization of J^,J^ and taking sum over all square bonds, we get (a -I- 
(3) (ri| S\n) =0 for periodic ladder and any ground state \ For the open ladder, one must add two boundary 
rung bonds with coefficients /? and a in order to obtain the last equation. So, if a 7^ — /3, then any ground state is a 
singlet. 

Consider now the ladder with site-dependent spin exchange and set again ai = ±f3i in (jA2|) . Note that if all 7/ 7^ 0, 
which makes strict the inequality sign in (jA5|) . then due to the discussions above the mean value of the spin of every 
rung vanishes for any ground state. Therefore, all ground states are singlets if Jj^ > \Ji^i\ + \Ji^\. In fact, the 
inequality is nonstrict for positive values of diagonal couplings. This becomes apparent if we perform the sum over 
square bonds instead of rung bonds using the parametrization ai = Pi- Then the ground-state mean value of the total 
spin of each such bond vanishes as has been mentioned above. Taking the sum over all square bonds for the periodic 
ladder and the sum over half of square bonds in checkerboard order for the open ladder with even number of rungs, 
we obtain the zero expectation for the total spin of the model provided that 

j,^>jr+jr-i, jr>o. (a?) 

For the open ladder with odd number of rungs, one of the boundary rung spins must be added to that sum. The 
ground state then is a singlet if the corresponding coefficient [71 or jpf in (|Aip ] is nonzero. This means that the 
inequality in one of the two boundary relations in (|A5[) must be strict, i.e., > or > J^_i- 
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